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1.  Introduction 


Inhibition  of  anion  exchange  reactions  during  molecular  beam  epitaxy  is  critical  to  the 
synthesis  of  mixed-anion  heterojunction  interfaces.  These  interfaces  are  the  basic  structural 
building  blocks  of  advanced  electronic  and  optoelectronic  devices  that  are  critical  to  DoD 
applications.  Anion  exchange  is  a  key  surface  and  subsurface  reaction  leading  to  interfacial 
intennixing.  Other  related  reactions  include  segregation,  diffusion,  and  clustering.  Interfacial 
intermixing  results  in  uncontrolled  strain  and  related  bandstructure  modifications  that  degrade 
device  performance,  reliability,  and  manufacturability. 

During  the  course  of  this  project,  we  have  studied,  both  experimentally  and  theoretically,  the 
fundamental  structural  properties  and  process  conditions  that  drive  anion  exchange  and  related 
reactions.  We  have  applied  our  techniques  to  As-,  P-,  and  Sb-based  mixed  anion 
heterostructures.  We  have  also  applied  advanced  modelling  techniques  to  enable  the  control  of 
anion  exchange  to  Sb/As  device  materials. 

1.1  Key  Findings 

The  key  finds  of  this  effort  are  listed  below. 

(1)  A  thermodynamic  approach  was  developed  to  assess  the  extent  of  anion  exchange  reactions 
during  the  heteroepitaxy  of  dissimilar  anion  III-V  compound  semiconductor  structures.  We 
found  that  the  extent  of  anion  exchange  can  be  predicted  by  the  change  in  the  Gibbs  free 
energy.  Bond  strength  changes  can  only  be  used  as  a  guide  in  comparing  the  relative  tendency 
for  exchange,  rather  than  as  a  criterion  for  exchange.  The  driving  force  for  anion  exchange 
strongly  depends  on  the  conditions  during  interface  formation.  A  number  of  important  factors, 
including  bond  strength,  misfit  strain  energy,  surface  structure  and  energy,  the  equilibrium 
between  dimers  V2  and  tetramers  V4,  and  segregation,  are  discussed  in  terms  of  their 
contributions  to  the  thermodynamics. 

(2)  A  wide  range  of  As/Sb,  Sb/As,  P/As  structures  exhibiting  different  degrees  of  anion 
exchange  reactions  were  synthesized  and  fully  characterized.  Superlattice  structures  were 
synthesized  by  allowing  a  dissimilar  anion  fluxes  to  impinge  on  a  binary  surface.  The 
intennixed  layers  are  characterized  using  high  resolution  x-ray  diffraction,  spectroscopic 
ellipsometry,  and  x-ray  photospectroscopy  (XPS).  The  importance  of  substrate  temperature, 
flux  exposure  time,  anion  surface  richness  (as  revealed  through  the  surface  reconstruction),  and 
incident  anion  molecular  species  was  examined.  Sb/As  structures  exhibit  Sb  surface 
segregation,  and  the  subsequent  incorporation  of  the  segregating  surface  population  into 
epitaxial  overlayers.  As/Sb  exchange  was  significant  and  enhanced  by  the  use  of  incident 
dimer  fluxes.  Competitive  reactions  between  exchange  and  isoelectronic  compound  formation 
yields  complex  structures  and  dependences  on  the  factors  examined.  P/As  exchange  was  also 
significant  and,  as  with  the  case  of  As/Sb,  competitive  reactions  were  found.  As-rich 
reconstructions  were  found  to  inhibit  exchange,  thereby  showing  the  importance  of  the  surface 
structure  to  subsequent  exchange.  Chemical  models  were  articulated  for  each  of  the  systems. 

(3)  A  hybrid  neural  network  model  was  constructed  by  characterizing  the  growth  of  GaAsi. 
yPy/GaAs  strained  layer  superlattices  (SLS)  grown  on  (001)  GaAs  substrates  by  molecular 
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beam  epitaxy.  These  heterostructures  were  formed  by  the  P2  exposure  of  an  As-stabilized 
GaAs  surface,  and  ex-situ  high-resolution  x-ray  diffraction  (HRXRD)  was  performed  to 
determine  the  phosphorus  composition  at  the  interfaces.  A  first-order  kinetic  model  was  then 
developed  to  describe  the  mechanisms  of  anion  exchange,  surface  desorption,  and  diffusion.  A 
semi-empirical  hybrid  neural  network  was  used  to  estimate  the  parameters  of  the  kinetic  model 
and  analyze  the  microscopic  processes  occurring  at  the  interfaces  of  the  mixed  anion  III-V 
heterostructures.  The  phosphorus  diffusion  process  in  GaAs  was  estimated  to  have  a  diffusion 
coefficient  of  D  =  7x  10"15  exp  (-0. 1  \/ksTs )  cmV1  for  samples  with  Pas4  =  4x  10'6  torr  and 
exhibited  enhanced  phosphorus  intennixing  for  samples  with  lower  As-stabilizing  fluxes. 

(4)  The  InAs-AlSb  high  electron  mobility  tTransistor  (HEMT)  is  of  critical  importance  for 
future  millimeter  wave  and  high-speed  electronic  devices.  We  have  achieved  record  300  K 
electron  mobility  in  InAs  channel  structures.  We  applied  various  growth  techniques  to  the 
critical  mixed-anion,  As/Sb,  inverted  interface  to  discover  the  impact  on  device-relevant 
properties. 

(5)  We  have  made  new  discoveries  regarding  the  clustering  and  dislocation  microstructure 

of  GalnP-GaAs  heterostructures  that  derive  from  anion  exchange  reactions.  The  infonnation  is 
contained  in  [9]  and  [10]  in  Section  4  and  not  included  in  the  body  of  this  report. 

(6)  In  addition,  while  outside  of  the  scope  of  our  statement  of  work  and  not  directly  supported 
by  this  effort,  we  applied  exchange  techniques  to  nanostructures  (quantum  dots)  to  show 
that  exchange  reactions  can  enhance  or  inhibit  quantum  dot  formation  and  can  lead  to 
improved  vertical  self-assembly. 
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2.  Detailed  Discussions 


2.1  Thermodynamic  Modelling 
2.1.1  Basic  Models 

Bond  strength  energies  are  often  used  to  predict  whether  anion  exchange  will  occur  or  not.  It  is 
generally  assumed  that  weaker-for-stronger-bond  exchange  (e.g.  Sb-for-As,  As-for-P,  and  Sb- 
for-P)  cannot  occur  under  typical  MBE  conditions  [1].  However,  this  was  proved  wrong  by 
Wang  [2],  who  reported  that  up  to  1  ML  of  As  can  be  replaced  by  Sb  when  the  As-tenninated 
InAs  (001)  surface  was  exposed  to  the  Sb2-flux  at  a  temperature  of  about  380  °C.  As-for-P 
exchange  was  observed  experimentally  much  earlier  [3].  A  detailed  study  showed  that  about 
two  monolayers  of  InAs  are  formed  on  the  InP  (001)  surface  as  a  result  of  As-for-P  exchange 
[4].  Further  studies  showed  that  As/P  exchange  can  also  occur  under  metal-organic  vapor 
phase  epitaxy  (MOVPE)  conditions  [5]. Our  results  in  Section  2.2  show  a  range  of  intermixing 
for  various  conditions,  as  well.  Apparently,  bond  strength  cannot  be  generally  used  as  the 
primary  criteria  for  predicting  whether  or  not  anion  exchange  will  occur. 

We  have  exploited  a  thennodynamic  approach  to  analyze  anion  exchange  reactions,  with  a 
focus  on  the  reaction  driving  force  and  direction  under  standard  MBE  conditions.  A  number  of 
important  factors,  including  bond  strength,  misfit  strain  energy,  surface  structure  and  energy, 
the  equilibrium  between  dimers  V2  and  tetramers  V4,  and  segregation,  were  discussed  in  tenns 
of  their  contributions  to  the  thermodynamics. 

Fig.  2.1.1  shows  a  simplified  model  of  anion  exchange  under  MBE  conditions.  The  surface  is 
exposed  to  a  VB2  flux  with  a  beam  equivalent  pressure  (BEP),  P0,  or  beam  flux,  Jo,  for  a  unit  of 
time.  Then,  the  VB2  beam  is  shut  off  and  the  system  is  relaxed  into  equilibrium.  After  the 
exchange  reaction,  VB2  molecules  depart  the  surface  with  a  flux,  JB2,  or  BEP,  PB2.  The 
exchange  product  for  the  group  V  species  may  be  a  mixture  of  Va2  and  Va4  (e.g.  AS2  and  AS4), 
depending  on  the  exchange  conditions. 


t 


t 


t 


Substrate:  IIIVA 


Sole  the  downward  arrow  represent  the  inward  flux,  and  the  upward  ones  the 
desorption  flux.  The  thK&ness  tebmiaDcall*  represents  the  magnitude  ol  the  fluxes 


Figure  2.1.1:  Simplified  thermodynamic  model  of  exchange  reaction. 


3 


Note  that  these  theoretical  conditions  for  the  thermodynamics  study  form  the  basis  for  our 
experimental  studies  in  Section  2.2  of  this  report.  For  clarity,  we  first  consider  Va2  only.  The 
flux  and  BEP  for  the  product  gas  Va2  are  Ja2  and  Pa2,  respectively.  During  exchange,  the 
surface  is  assumed  to  remain  flat,  smooth,  and  coherent  to  the  substrate.  This  is  a  simplifying 
assumption  that  is  not  generally  realized  as  the  lattice  strain  drives  morphological  transitions. 
Bulk  diffusion  is  neglected,  but  in  the  exchanged  layer,  diffusion  is  presumed  to  be  fast  enough 
not  to  interfere  with  the  exchange  process.  Assuming  no  deposition  or  sublimation,  and 
neglecting  the  background  pressures,  mass  conservation  requires  that 


Jo  -  Jb2  -  Ja2 


(1) 


The  anion  exchange  is  described  as  an  interfacial  reaction, 


III-Va(s)  +  1/2  VB2(g)  ->  HI-Vb(s)  +  l/2VA2(g)  (2) 


The  ‘s’  and  ‘g’  in  parentheses  stand  for  solid  and  gas,  respectively.  The  equilibrium  constant, 
K,  for  this  reaction  can  be  shown  to  be 


K  -  Kiii-va/Kih-vb 


(3) 


where  Kih_va  and  Khi-vb  are  equilibrium  constants  for  sublimation 


III-V(s)  ^  111(g)  +  l/2V2(g) 


(4) 


Listed  in  Table  2.1  are  the  equilibrium  constants  for  the  III-V  compounds,  which  are  calculated 
based  on  thermodynamic  data  [6],  except  those  for  As2.  As  pointed  by  Hurle  [7],  the  standard 
heat  of  formation  of  As2  molecules  (AH°298  ~  53.5  kcal/mole)  is  too  large.  We,  therefore,  use 
the  data  provided  by  Pupp  [8],  For  As2,  the  fonnation  heat  is  44.4  kcal/mol,  and  34.4  kcal/mol 
for  AS4.  The  equilibrium  constant  for  As2  <-»  AS4  is  also  included  in  Table  2.1.  We  will  use  this 
value  in  the  discussion  of  the  effect  of  AS4  on  the  exchange  reaction. 


4 


Table  2.1:  Equilibrium  constants  for  sublimation  of  III-V  Compounds  and  the  maximum 
sublimation  temperature  (TSUbi)  and  the  temperature  of  non-congruent  dissociation 

(Ts„b2) 


Reaction 

(KT  in 

Taohl  CO 

1*ib2  C°C) 

Aiaj+i/iPai  -Alt',* 

3.86x10'  -exp(-S,S4/kT) 

Alfg)  t  l/2AS2(gj  *-*■  AiASja) 

Alja+l/SSb^g)  -A1SV) 

5,00X1 0lIexK-S.58TO 

2 .7  9x  1 G1  -  ex.p(  -  S.  1 0/kT) 

W2 

974 

Ga.£.+].J|2P''isi  *  *  GaPiE) 

4.19xlOnexp(“4.5SitT) 

571 

774 

Ga^+l/IAs^  *-*■  GaAsw 

Ga^+ 1  /ISb^gj  —  GaSbfij 

5  ,44x  1 Q1  Lexp(  -L55;kT> 
2,71x1 0i:'exp(  — 4.53/fcI) 

630 

723 

-  XnP(>;i 

l.bbxlO^expC  4.14,lcT) 

268 

634 

lut  1  /2 A  sag]  ~  liiAsts) 
la^+ljCSb^.  -  InSb,:t;. 
Ga^-l-l/4SMs]  "  GaSbitj 
*+  lnSb(,) 
As^eI  w  3/2As4^ 

2 .05  X 1  e?ip(  -  4,09/kT) 

4 ,37x  1 01Qexp(  -3.97;kTj 
2,73x1 01  °exp(  -  3 . SS/kT) 
4.39x1  tfejcj!(---3i32/lT) 

5 ,00x  1 0“  sejrp(  1 . 1 8/kT) 

508 

m 

The  product  III-VB  mixes  with  III-VA  and  forms  a  III-VA(i-x)Vbx  solid  solution.  This  process  is 
spontaneous  for  all  the  III-V  compound  semiconductor  alloys  in  the  MBE  temperature  range, 
although  the  solution  may  be  metastable,  as  is  the  case  for  Ga-As(i.x)Sbx  ternary  alloys.  Based 
on  a  regular  solution  model  [9],  the  activities  of  the  components  in  the  solution,  denoted  by 
[III-V]  take  the  form: 


[III-VA]  =  x  exp[Q(l  -  x)2/kT]  (5) 

[III-Vb]  =  ( 1  -  x)  exp[Qx2/kT]  (6) 


where  the  interaction  energy  Q  and  kT  are  measured  in  electron  volts  (eV). 

Let  ps  represent  the  atomic  density  on  the  surface  and  N  represent  the  number  of  monolayers 
participating  the  exchange  reaction,  the  total  number,  Nt,  of  VA  atoms  being  replaced  by  VB  is 

Nt=  xNps  =  xNs  =  2Ja2  (7) 


By  combining  Eqs.  (1),  (7),  we  get 
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where  mA2  and  mB2  are  the  molecular  weights  of  VA2  and  Vb2-  When  the  reaction  reaches  its 
equilibrium,  Q  =  K,  the  molar  fraction  of  exchange,  x,  is  determined.  In  this  equation,  the 
exponential  term  measures  how  far  the  solid  solution  deviates  from  an  ideal  solution. 

The  number  of  monolayers,  N,  in  the  exchange  reaction  is  unknown.  In  our  calculation,  we  use 
the  ratio  Jo/Ns  as  an  input.  One  can  readily  incorporate  the  thickness  of  the  exchanged  layer 
into  our  results.  However,  for  the  sake  of  clarity  and  simplicity  (to  compare  different 
reactions),  our  discussion  assumes  one  monolayer  exchange  and  the  exposed  surface  is  (001). 
In  this  case,  for  BEP  =  10'6  Torr,  Jo/Ns  is  about  0.2  (ps~5.5  x  1014/cm2). 

We  consider  P-for-As  (strong-for-weak)  and  Sb-for-As  (weak-for-strong)  as  examples  to 
illustrate  the  thermodynamic  characteristics  of  the  exchange  reactions.  As  will  be  shown 
below,  the  equilibrium  constants  for  the  exchange  reactions  are  weakly  temperature-dependent. 
The  impact  of  temperature  is  much  less  remarkable  than  other  factors,  such  as  J0.  In  our 
examples,  we  have  tried  to  use  the  same  temperatures  as  experimentally  used. 

2.1.2.  P/As  Exchange 

Strong-for-weak  bond  exchange  is  expected  to  occur  under  normal  exposure  experimental 
conditions.  This  is  what  we  obtain  from  thermodynamic  prediction,  as  shown  in  Fig.  2.1.2  (a) 
and  (b),  where  the  exposed  surfaces  are  GaAs  (001)  and  InAs  (001),  respectively. 
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T  =  800K 


Figure  2.1.2:  Exposure  to  P2.  (a)  GaAs  (001)  exposed  to  P2  and  (b)  InAs  (001)  exposed  to 

P2. 


The  arrowheads  indicate  the  exchanged  molar  fraction  at  equilibrium.  The  exchange  reactions 
and  their  equilibrium  constants  for  these  two  cases  are: 

GaAs(s)  +  1/2  P2(g)  ->  GaP(s)  +  l/2As2(g),  K  =  1 .30  exp(0.035/kT)  (9) 


InAs(s)  +  1/2  P2(g)  ->  InP(s)  +  l/2As2(g),  K  =  1.23  exp(0.052/kT)  (10) 


It  can  be  seen  that  the  exchanged  molar  fraction,  x,  of  GaP  or  InP  increases  with  the  ratio  of 
J0/Ns.  Under  small  fluxes  (J0/Ns  <  ~  0.25),  almost  all  the  incoming  P2  molecules  are  exchanged. 
When  the  flux  increases  to  0.5NS,  which  corresponds  to  a  BEP  of  about  1.6  x  10'6  Torr 
(assuming  one  monolayer  exchange),  about  60%  of  the  As  atoms  in  an  InAs  (001)  monolayer 
are  replaced  by  P  atoms.  It  can  also  be  seen  that  the  degree  of  exchange  increases  with  the 
increase  in  bond  strength  difference. 
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It  should  be  noted  that  the  backward  reaction  represents  As-for-P  (weak-for-strong)  exchange. 
Therefore,  under  the  same  J0/Ns  condition,  the  As-for-P  exchange  can  hardly  occur  for  small 
J0/Ns.  To  get  appreciable  exchange,  the  J0/Ns  should  be  larger  than  ~  1  x  10"6  Torr. 

2.1.3.  Sb/As  Exchange 

In  Fig.  2.1.3,  the  variation  of  the  reaction  quotients  with  molar  fraction,  x,  of  exchange  under 
different  J0/Ns  conditions  is  shown. 


Figure  2.1.3:  Sb2  Exposure,  (a)  GaAs  (001)  exposed  to  Sb2  and  (b)  InAs  (001)  exposed  to 

Sb2. 


The  exposed  surfaces  are  GaAs  (001)  and  InAs  (001).  The  exchange  reactions  and  equilibrium 
constants,  respectively,  are: 


GaAs(s)  +  1/2  Sb2(g)  ->  GaSb(s)  +  l/2As2(g),  K  =  0.20  exp(-0.02/kT)  (11) 
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InAs(s)  +  1/2  Sb2(g)  ->  InSb(s)  +  l/2As2(g),  K  =  4.69  exp(-0. 12/kT) 


(12) 


As  expected,  the  weak-for-strong-bond  exchange  is  very  small.  For  example,  under  J0/Ns  = 
0.25,  which  corresponds  to  a  BEP  of  1.6  x  10'6  Torr,  x  =  0.05  for  GaAs  (001).  In  contrast, 
almost  all  As2  molecules  are  exchanged  for  the  As-for-Sb  reactions  under  the  same  conditions, 
qualitatively  in  agreement  with  the  experimental  observations. 

In  reality,  there  are  a  number  of  factors  that  may  have  large  impact  on  the  thermodynamic 
equilibrium.  They  may  not  be  able  to  reverse  the  reaction  directions,  but  they  are  indeed  able  to 
shift  the  equilibrium  point.  These  factors  are  discussed  in  detail  below. 

2.1.4.  Bond  Strength 

Although  the  exchange  reaction  is  often  discussed  in  terms  of  bond  strength,  the  direction  of 
exchange  reaction  is  detennined  by  the  total  free  energy  change,  rather  than  by  bond  strength. 
However,  under  normal  MBE  conditions,  bond  strength  can  be  used  as  a  first  guide  for 
determining  the  tendency  of  exchange.  Fig. 2.4  shows  the  equilibrium  constants  in  the  MBE 
temperature  range  for  some  typical  exchange  reactions.  Note  the  weak  temperature  dependence 
of  the  equilibrium  constants.  It  can  be  seen  that  the  larger  differences  in  bond  strength  have 
larger  equilibrium  constants  for  strong-for-weak-bond  exchange  (or  smaller  K  for  weak-for- 
strong-bond  exchange),  implying  larger  degree  of  exchange. 


Figure  2.4:  Equilibrium  constant.  Temperature  dependence  of  equilibrium  constant  for 

some  typical  exchange  reactions 


2.1.5.  Mixing  Interaction  Energy 

When  foreign  or  solute  atoms  are  put  into  a  host  lattice,  lattice  distortion  is  produced 
corresponding  to  the  atomic  size  difference.  The  amount  of  strain  energy,  related  to  the  lattice 
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distortion,  and  also  the  chemical  potential  of  the  constituent  components,  is  measured  in  tenns 
of  the  mixing  interaction  energy.  The  effect  of  the  interaction  on  the  equilibrium  of  the 
exchange  reaction  is  reflected  by  the  term  exp[(l-2x)  Q/kT)].  The  interaction  energy,  Q,  which 
is  related  to  the  regular  solution  model,  has  been  shown  to  be  directly  proportional  to  the 
square  of  lattice  mismatch.  Based  on  the  delta-lattice-parameter  model  proposed  by 
Stringfellow  [10],  Cl  is  147.5  meV  for  Ga-As-Sb,  98.5  meV  for  In-As-Sb,  42.5  meV  for  Ga- 
As-P,  and  25.2  meV  for  the  In-As-P  system.  Consequently,  the  greatest  effect  of  the  interaction 
energy  is  expected  in  the  Ga-As-Sb  ternary  alloys.  The  equilibrium  point  is  shifted  in  the  lower 
x  direction  when  x  <  0.5  (usually  for  weak-for-strong-bond  exchange),  and  for  x  >  0.5,  to  the 
higher  x  direction.  Interpretation  of  this  effect  is  simple.  When  x  <  0.5,  the  lattice  distortion 
energy  increases  with  addition  of  solute  Vb  atoms  to  the  host  III-Va  lattice.  This  energy 
increase  process  automatically  results  in  shifting  of  the  equilibrium  point  to  smaller  x.  For  x  > 
0.5,  the  host  lattice  is  converted  from  III-Va  to  III-Vb.  Increasing  x  is  equivalent  to  decreasing 
the  solute  (now  Va)  concentration,  and  hence  the  equilibrium  is  shifted  to  higher  x  values.  A 
quantitative  comparison  is  shown  in  Fig.  2.1.5.  It  can  be  seen  that  remarkable  effects  can  be 
expected  when  the  equilibrium  x  is  away  from  0.5. 


Molar  fraction  of  exchange,  x 


Figure  2.1.5:  Surface  Exposure  to  Sb2  and  P2.  (a)  InAs  (001)  exposed  to  Sb2;  (b)  GaAs 

(001)  exposed  to  P2. 


2.1.6.  Strain  Energy 

As  a  result  of  atomic  differences,  the  lattice  constant  of  the  solid  solution  III-VA(i-x)Vbx  is 
different  from  that  of  the  substrate  III-Va,  giving  rise  to  lattice  mismatch.  Consequently,  the 
surface  layer  is  strained  and  elastic  strain  energy  is  built  up  during  the  anion  exchange,  thereby 
increasing  the  chemical  potential  of  the  constituent  components.  The  amount  of  strain  energy 
stored  depends  on  the  molar  fraction,  x.  Assuming  the  lattice  constants  obey  Vegard’s  law  and 
no  plastic  relaxation  occurs,  the  strain  energy  Estr  (in  eV)  takes  the  fonn  of 
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where  tn,  is  atomic  weight,  p  the  density  (g/cm  );  the  C’s  are  elastic  constants  (Pa),  s0  the 
misfit  strain  when  x  =  1 . 

As  expected,  the  stored  strain  energy  in  the  exchanged  surface  layer  always  shifts  the 
equilibrium  point  to  lower  molar  fractions.  In  other  words,  the  strain  effect  always  resists  the 
exchange  process  on  a  strain- free  flat  surface.  This  is  in  contrast  with  the  case  where  anion 
exchange  reduces  strain  for  an  initial  three-dimensional  (or  islanded)  surface.  In  this  case,  the 
strain  energy  in  the  system  maybe  reduced  by  the  anion  exchange,  thereby  providing  a  driving 
force  for  size  reduction  of  three-dimensional  islands,  or  flattening  of  the  surface. 

The  effect  of  strain  energy  on  the  equilibrium  is  shown  in  Fig.  5,  where  the  contributions  of 
non-ideality  of  the  solid  solution  and  strain  effect  are  compared  with  respect  to  the  ideal 
solution.  More  comparisons  are  shown  in  Figs.  2.1.2,  2.1.3,  and  2.1.6. 


Figure  2.1.6:  InAs  (001)  exposed  to  Sb2. 


For  group-V  species,  an  equilibrium,  exist  between  tetrameric  and  dimeric  form: 

v2(g)  <->  i/2V4(g),  kV4  =  pf;p;l  (14) 


Depending  on  the  exchange  temperature,  this  equilibrium  can  have  drastic  impact  on  the 
exchange  reaction.  Within  the  temperature  range  of  530-730  °C,  P2  is  the  dominant  phosphorus 
species  over  InP.  Only  up  to  10%  of  the  total  phosphorus  flux  consists  of  P  tetramers  [1 1,12], 
Fluxes  originating  from  the  (001)  InP  surface  heated  in  the  UHV  between  330-420°C  consist 
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of  only  P  dimers  [11].  These  results  indicate  that  in  the  practical  exchange  temperatures,  the 
involvement  of  P  tetramers  can  be  neglected. 


For  Sb2(g)  <-»  l/2Sb4(g),  it  is  predicted  that  the  dominant  equilibrium  antimony  species  is  the 
Sb  tetramer  below  530°C,  based  on  the  thermodynamic  data  [13].  Under  a  total  pressure  in  the 
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range  of  10'  -  5x10'  Torr,  about  90%  of  the  molecules  are  Sb4.  Therefore,  Sb4  molecules  are 
expected  to  dominate  in  the  flux  from  the  surface  at  typical  MBE  conditions. 

For  As2(g)  <-»  l/2As4(g),  experimental  results  indicate  that  below  350°C,  the  tetrameric  As4 
dominates  in  the  flux  from  a  heated  GaAs  surface,  irrespective  of  the  type  of  arsenic  molecules 
impinging  on  the  surface.  At  temperatures  of  450°C  and  higher,  dimeric  arsenic  molecules  take 
over  in  dominance  [14,15].  These  results  essentially  agree  with  the  thermodynamic  predictions 
[16,17]. 
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Fig.  2.1.6  shows  the  effect  of  As4  involvement  in  the  exchange  reaction  when  InAs  is  exposed 
to  Sb2.  For  convenience  of  comparison,  the  no-As4-involement  data  is  also  duplicated  form 
Fig.  2.1.3a.  It  can  be  seen  that  the  effect  of  As4  on  the  equilibrium  is  remarkable.  The  molar 
fraction  x  increases  from  0.30  to  0.66  under  the  same  exposing  conditions.  From  Fig.  2.1.4,  it 
can  also  been  seen  that  the  equilibrium  constant  is  very  strongly  temperature  dependent. 


2.1.7.  Surface  Structure  and  Surface  energy 

Of  all  the  factors  influencing  the  reaction  thermodynamics  and  kinetics,  surface  structure  and 
surface  energy  may  be  the  most  important.  Semiconductors  have  reconstructed  surface 
structures  in  order  to  assume  a  lower  energy  state  under  MBE  conditions.  Therefore,  the 
surfaces  may  have  considerably  different  thermodynamic  properties  from  their  bulk  state.  As  a 
result,  the  equilibrium  state  that  is  predicted  based  on  the  bulk  material  properties  may  be  far 
away  from  the  actual  one.  Unfortunately,  few  thermodynamic  data  are  available  for  any 
quantitative  discussion.  Intuitively,  the  shifting  of  the  equilibrium  point  is  expected  to  be  more 
significant  for  the  case  where  anion  exchange  results  in  a  transition  in  surface  reconstruction 
(e.  g.  (1  x  3)  GaSb  to  (2  x  4)  GaAs)  than  the  case  where  no  surface  structure  transition  is 
involved. 
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Anion  exchange  mostly  occurs  in  the  first  one  or  two  monolayers.  Other  reactions  lead  to  much 
more  broadened  interfaces.  Therefore,  the  surface  energy  occupies  a  large  portion  in  the  total 
energy  involved  in  the  exchange  process.  Even  a  small  change  in  the  surface  energy  may  give 
rise  to  a  large  shift  in  the  equilibrium  state.  However,  due  to  the  lack  of  surface  energy  data,  we 
can  only  make  some  qualitative  assertions.  Consider  the  As-terminated  (001)  surface  energy.  It 
is  reported  that  this  surface  energy  is  2100  meV  per  mrT  [18].  So  for  each  As  atom,  surface 
energy  is  about  335  meV,  up  to  ten  times  higher  than  the  misfit  strain  energy.  For  weak-for- 
strong-bond  exchange  (e.g.  Sb-for-As),  the  surface  energy  will  be  reduced  because  of  the 
formation  of  weak  bonds  on  the  surface,  and  hence,  the  exchange  reaction  will  be  favored 
energetically.  In  contrast,  strong-for-weak-bond  exchange  will  tend  to  increase  the  surface 
energy,  and  hence,  will  be  reduced  by  the  change  in  surface  energy. 


2.1.8.  Segregation 

In  our  model,  the  exchanged  molecules  are  presumed  to  depart  the  surface  after  exchange.  This 
is  not  always  true.  Segregation  can  occur  in  some  material  systems.  Although  this  is  a  post¬ 
exchange  process,  it  can  have  large  impact  on  the  equilibrium  point.  The  most  problematic 
anion  is  antimony  (Sb).  It  has  been  reported  that,  in  the  heteroepitaxy  of  InAs/GalnSb 
structures,  the  InAs-on-GalnSb  interfaces  are  tens  of  monolayers  thick.  The  presence  of  Sb  on 
the  growing  surface  can  be  detected  even  at  an  InAs  layer  thickness  of  200  mn  [19].  Much 
more  convincing  results  are  reported  recently  by  Kaspi  [20],  who  used  desorption  mass 
spectrometry  to  monitor  the  chemistry  of  the  flux  departing  the  surface  in-situ.  The  segregated 
Sb  is  incorporated  into  the  lattice  of  the  InAs,  forming  a  broad,  strained  and  Sb-graded 
interface. 

2.2.  Chemical  Models  for  Sb/As,  As/Sb,  and  P/As  Heterostructures 

We  show  in  this  section  that  their  significant  complexity  in  the  anion  exchange  and 
interdiffusion  process.  In  comparison  to  the  previous  section,  we  see  that  surface  structure  and 
anion  interactions  are  key  to  our  experimental  observations.  Revisiting  of  the  analysis  in 
Section  4. 1  is  advisable  from  these  results. 

The  exposure  of  GaAs  and  GaSb  surfaces  to  P,  As,  and  Sb  beams  was  carried  out  to  understand 
the  extent  of  and  driving  forces  for  anion  exchange  during  structure  fonnation.  Numerous 
factors  were  compared.  These  include:  substrate  temperature,  surface  preparation 
(reconstruction),  anion  flux  species  (dimer  or  tetramer)  and  exposure  time.  Superlattice 
structures  are  prepared  not  through  direct  growth,  but  by  the  exposure  of  the  semiconductor 
surface  to  a  different  anion  flux  at  a  specific  temperature  and  for  a  specific  time  (same  type  of 
structure  analyzed  theoretically  in  Section  2.1  of  this  report).  Through  the  analysis  of  data 
derived  from  a  range  of  characterization  techniques,  including  RHEED,  high  resolution  x-ray 
diffraction  (XRD),  spectroscopic  ellipsometry  (SE),  and  AFM,  we  have  created  models  of  the 
structures  and  inferred  primary  driving  forces  for  creating  the  observed  structures.  Chemical 
models  are  then  articulated  which  capture  the  reactions  present  during  junction  synthesis. 

In  our  studies,  As/Sb  and  Sb/As  anion  exchange  reactions  were  investigated  by  exposure  of 
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GaSb  surfaces  to  AS2  and  AS4  species,  and  by  exposure  of  GaAs  to  Sb2,  respectively.  It  is 
found  that  when  GaSb  surfaces  are  exposed  to  As,  competitive  reactions  occur.  These  reactions 
are  complex  and  indicate  why  a  thennodynamic  analysis,  while  important  as  “setting  the  stage” 
for  understanding  exchange,  do  not  accurately  predict  experimentally-observed  heterojunction 
characteristics.  As  shown  herein,  As  exchanges  with  Sb  in  the  GaSb  resulting  in  the  formation 
of  GaAs  bonds,  and  also  reacts  with  Sb  to  form  AsSby  isoelectronic  compounds.  The  relative 
degrees  to  which  GaAs  and  AsSby  are  formed  depend  on  the  temperature,  As  exposure  time, 
and  the  As  species  (i.e.,  AS2  or  AS4)  incident  on  the  surface.  It  is  found  that  AS2  is  more  likely 
to  exchange  with  Sb  in  comparison  to  AS4,  inhibiting  the  formation  of  AsSby.  AsSby  interface 
segregation  inhibits  continued  exchange.  A  specific  AS4  exposure  time  and  temperature  exists 
wherein  the  formation  of  AsSby  is  minimized.  AS4  and  AS2  have  been  previously  observed  to 
react  differently  on  GaSb  surfaces  [19].  We  also  observe  that  AS4  is  less  reactive  than  AS2,  and 
consequently  As  incorporation  into  GaSb  layers  is  reduced  when  As4  impinges  on  the  GaSb 
surface.  In  the  case  of  the  Sb2/GaAs  system,  Sb  for  As  anion  exchange  does  not  occur  to  a 
significant  degree,  and  surface  segregation  of  Sb  is  found,  with  the  subsequent  incorporation  of 
segregated  Sb. 


2.2.1  P/As:  Impact  of  Temperature,  Exposure  Time,  and  Surface 
Reconstruction 

As/P  intermixing  has  also  been  studied  for  some  time,  particularly  during  the  early  exploration 
of  the  synthesis  of  GalnAs/InP  by  chemical  beam  epitaxy.  Exposure  of  As-based  compound 
surfaces  to  P  is  a  standard  procedure  in  beam  switching  sequences  for  the  growth  of  a  number 
of  important  hetero structures.  When  GaAs  is  exposed  to  P,  the  group  V  element  may  undergo 
P  for  As  anion  exchange,  with  the  formation  of  a  thin  layer  of  GaP  or  GaPyAsi_y.  Degradation 
of  the  interface  quality  has  been  found  as  a  result  of  the  interfacial  segregation  and  exchange 
reaction  between  As  and  P.  [21] 

Both  (2x1  O'6  T)  or  high  (4x1 0‘6  T)  As4  beam  equivalent  pressures  were  using  during  the 
synthesis  of  GaAs  to  prepare  both  As-lean  and  As-rich  (relatively)  GaAs  surface 
reconstructions.  It  is  found  that  the  As-rich  reconstruction  reduces  the  phosphorus  surface 
coverage  and  hinders  P-for-As  exchange.  The  effect  of  surface  reconstruction  on  exchange  and 
the  functional  dependences  of  P  incorporation  on  temperature  and  exposure  time  indicate  that  P 
in-diffusion  does  not  control  the  anion  exchange.  We  introduce  a  reactive  process  that  involves 
As  on  the  GaAs  surface.  This  process  is  the  competitive  formation  of  AsP  compounds.  In  this 
framework,  As  on  the  surface  is  a  sink  for  P,  and  thereby  controls  the  extent  of  anion 
exchange.  A  chemical  model  including  this  competitive  reaction  and  AsP  isoelectronic  by¬ 
product  formation,  in  addition  to  P  diffusion  into  the  GaAs,  is  proposed.  It  can  be  seen  that  the 
competition  between  isoelectronic  compound  formation  and  exchange  is  a  generalizable 
behavior  across  the  materials  systems  we  have  investigated. 

Twenty-period  GaAsySbi_y/GaSb  (exposure  of  GaSb  to  As)  and  GaSbyAsi_y/GaAs  (exposure  of 
GaAs  to  Sb)  structures  were  grown  in  a  Varian  Gen  II  MBE  system  equipped  with  valved 
cracker  As  and  Sb  sources.  For  the  experiments  in  which  GaSb  was  exposed  to  As,  the 
structure  was  initiated  with  a  GaSb  buffer  layer  (-270 A)  deposited  onto  a  GaSb  substrate.  AS2 
and  AS4  fluxes  were  delivered  to  the  Sb-stabilized  GaSb  surface  for  10-30  s,  over  a  substrate 
temperature  range  of  420  °C  -  470  °C.  The  beam  equivalent  pressures  (BEP)  were  As4=4xl0’6 


14 


Torr,  of  AS2  =2x1  O’7  Torr  (the  BEP  changes  were  a  consequence  of  the  variation  in  the  cracker 
tip  temperature,  all  other  conditions  were  the  same). 

Ten-period  GaAsyPi.y/GaAs  (GaAs  exposure  to  P)  were  synthesized  in  a  Riber  2300  solid 
source  MBE  system  with  P  and  As  valved  cracker  sources.  During  GaAs  growth,  two  different 
As  BEPs,  4x1  O'7  T  and  2x1 0~6  T,  were  used  to  prepare  different  GaAs  surfaces.  RHEED 
measurements  showed  that  the  higher  As  flux  yielded  different  variants  of  the  (2x4)  As-rich 
surface  reconstruction,  while  the  lower  As  flux  results  in  an  (2x1)  GaAs  surface  reconstruction. 
The  thickness  of  the  GaAs  layer  was  approximately  200A.  P  incorporation  in  the  structure  was 
achieved  by  P2  exposure  (P2  BEP  =  4.5xl0‘6  Torr)  of  the  GaAs  surface  and  subsequent 
overgrowth  of  GaAs.  The  substrate  temperature  was  varied  from  420°C  to  520°C,  and  the  P2 
exposure  time  was  varied  from  10s  to  50s. 

For  the  Sb/As  samples,  angle-resolved  XPS  analysis  was  perfonned  at  take-off  angles  of  20°, 
45°  and  90°,  corresponding  to  sampling  depths  of  approximately  29,  61  and  87  A,  respectively. 
High-resolution  spectra  (pass  energy  =  17.90  eV)  of  the  Ga  3d,  As  3d,  Sb  2d  and  Sb  3d 
photoelectron  core  levels  were  acquired.  A  spin-orbit  splitting  separation  energy  of  9.34  eV 
exists  between  the  Sb3d3/2  and  Sb3ds/2  levels.  However,  the  Sb  3ds/2  peak  suffers  from 
interference  with  the  Ols  signal  at  the  binding  energy  (BE)  of  531  eV;  hence,  the  Sb  2d  core 
level  was  also  considered,  although  the  low  splitting  energy  of  1.2  eV  for  the  Sb2d3/2  and 
Sb2d5/2  renders  the  fit  analysis  difficult.  In  addition,  care  was  taken  since  the  Ga  3d  plasmon 
loss  tail  interferes  with  the  Sb  2d  peak.  Therefore,  bond  configurations  of  Sb  were  established 
by  analysis  of  both  Sb  3ds/2  and  the  Sb2d  signals.  Charging  effects  were  corrected  using  the  C 
Is  photoelectron  peak  from  surface  contaminants  at  285.0  eV.  The  near-surface  composition 
was  estimated  from  the  peak  areas  of  the  high-resolution  spectra  and  using  atomic  sensitivity 
factors  provided  by  Perkin-Elmer.  For  the  P/As  samples,  high  -resolution  spectra  of  the  Ga  3d, 
As  3d,  P  2p  O  Is  and  C  Is  photoelectron  core  levels  were  acquired  with  a  pass  energy  of  19.75 
eV.  The  position  of  all  the  above  peaks  was  calibrated  using  C  Is  with  binding  energy  BE  = 
284.6  eV.  Angle -resolved  XPS  measurements  were  perfonned  as  described  for  the  Sb/As 
samples. 

The  layered  structure  pseudodielectric  functions,  <s>  =  <Si>  +  i<S2>,  were  measured  in  the  1.5 
-5.5  eV  photon  energy  range  with  an  energy  resolution  of  0.02  eV  with  a  phase  modulated 
spectroscopic  ellipsometer  (UVISEL-Jobin  Yvon)  at  multiple  angles  of  incidence  in  the  range 
55°  -  80°  for  Sb/As  structures.  Correspondingly,  for  the  P/As  structures  the  0.75-6  eV  photon 
energy  region  with  an  interval  of  0.02  eV  was  used.  The  spectra  were  analyzed  using  the 
Bruggemann  effective  medium  approximation  [22].  The  experimental  data  were  fitted  by 
minimizing  the  squared  differences  between  the  measured  and  calculated  ellipsometric 
quantities  [23].  Layer  thicknesses  and  constituent  volume  fractions  were  used  as  the  fit 
variables.  Layer  compositions  are  evaluated  through  analysis  of  the  amplitude  and  energy 
position  of  the  Eo  and  Ei  interband  critical  points  (CPs).  To  obtain  a  quantitative  and  consistent 
analysis,  modeling  of  both  XRD  and  SE  spectra  (for  the  P/As  structures  only)  has  been 
performed  employing  a  multilayer  model  of  the  epitaxial  structure.  Specifically,  the  measured 
pseudodielectric  functions  of  the  samples  have  been  analyzed  using  a  repetition  (superlattice) 
model  GaAs-substrate/(GaAs/GaPyAsi_y)io/GaAs-oxide  (where  10  is  the  number  of  periods 
grown)  to  take  into  account  the  effect  of  the  surface  oxide.  A  linear  grade  approximation  was 
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used  to  model  the  in-diffusion  of  phosphorus.  The  linear  gradient  provides  good  fits  to  both  SE 
and  XRD  spectra,  and  it  was  verified  that  the  use  of  an  exponential  diffusion  gradient  does  not 
improve  the  goodness  of  fit  of  the  XRD  data. 

The  structures  were  also  characterized  by  high-resolution  X-ray  diffraction  (XRD)  in  order  to 
derive  their  quality  and  structural  parameters  such  as  the  thickness  and  composition  of  layers. 
XRD  analysis  was  performed  with  a  Philips  MRD  system  using  Cu  Ka  radiation  to  obtain 
(002)  9/29  scans.  The  superlattice  structures  were  simulated  with  the  dynamical  diffraction 
theory  using  the  Bede  RADS  software. 

Finally,  all  of  the  P/As  and  a  selection  of  the  As/Sb  and  Sb/As  structures  were  measured  using 
atomic  force  microscopy  (AFM)  and  all  were  monitored  during  growth  by  RHEED. 

2.2.2.  Sb/As  and  As/Sb  Reactions 

Figure  2.2.1  compares  the  Ga  3d,  Sb  2d,  Sb  3d5/2  and  As  3d  photoelectron  peaks  acquired  for 
GaSb  and  GaAs  substrates,  and  for  a  GaSb  surface  exposed  to  As2  and  AS4  fluxes  for  19s  at 
420  °C.  In  the  Ga  3d  signal  of  the  reacted  samples,  the  GaAs  (BE  =  19.2  eV)  and  GaSb  (BE  = 
19.6  eV)  fit  components  have  not  been  resolved  because  of  the  close  binding  energy  values, 
and  have  been  included  in  the  same  peak  at  19.5  eV,  which,  therefore,  has  a  relatively  large  full 
width  at  half  maximum  (FWHM)  of  1.3  eV.  Because  the  analysis  is  ex-situ,  the  Ga203  fit 
component  at  20.5  eV  is  also  present.  In  the  Sd  4d  photoelectron  spectra,  the  Sb203  and  Sb205 
contributions  have  also  been  included  in  the  same  fit  because  of  the  small  difference  in  binding 
energies.  The  splitting  energy  between  Sb  2d3/2  and  Sb2ds/2  has  been  fixed  at  1 .2  eV,  while  the 
energy  separation  between  Sb  in  GaSb  and  Sb  in  Sb203  and  Sb205  has  been  set  at  2.5  ±  0.4  eV. 
In  the  As  3d  spectra  of  the  GaSb  samples  exposed  to  the  arsenic  flux,  a  high  energy  component 
appears  in  addition  to  that  of  GaAs  and  As-oxides.  This  high  energy  peak  has  been  ascribed  to 
the  formation  of  isoelectronic  compound,  AsSby,  during  the  As  for  Sb  anion  exchange  reaction. 
The  AsSby  fit  component  is  very  broad  because  it  can  include  isoelectronic  components  of 
different  stoichiometry.  Specifically,  As  can  react  with  Sb  resulting  in  AsSb,  As2Sb3,  As3Sb. 
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Figure  2.2.1:  Photoelectron  Peaks.  Comparison  of  Ga3d,  Sb2d,  Sb3d5/2  and  As4d  high 
resolution  photoelectron  peaks  for  GaSb  and  GaAs  substrates  and  for  two  typical  GaSb 
layers  exposed  to  AS2  at  10s  and  to  AS4  at  30s  and  420°C. 


The  high  binding  energy  of  AsSby  in  the  As3d  implies  that,  in  the  AsSby  compound,  the  Sb  is 
electronegative,  and  the  As  is  electropositive.  Consequently,  the  AsSby  contribution  should  be 
at  low  binding  energy  in  the  Sb2d  photoelectron  spectrum,  and  its  binding  energy  value  could 
be  quite  close  to  that  of  GaSb.  This  would  be  consistent  with  the  increase  of  the  ratio  of  the 
low  binding  energy  peaks  to  the  Sb  oxide  peak  for  samples  with  a  corresponding  high 
contribution  of  the  AsSby  in  the  As3d  peak.  The  angle  resolved  analysis  of  the  GaSb  layer 
exposed  to  the  AS4  flux  for  10s  at  420°C  is  shown  in  Fig.  2.2.2,  and  shows  that  the  intensity  of 
the  AsSby  component  increases  with  take-off  angle,  i.e.,  the  sampling  depth,  and  become 
comparable  to  the  GaAs  component  at  high  take-off  angle,  indicating  that  AsSby  segregates  at 
the  GaAsySbi-y/GaSb  interface. 
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Figure  2.2.2:  Angle-resolved  As3d  photoelectron  peak  recorded  for  a  typical  GaSb  layer 
exposed  to  AS4  for  10s  at  420°C. 


The  formation  of  GaAs  bonds  as  a  product  of  exchange  and  AsSby  formation  depends  on  the 
incident  As  species,  i.e.,  AS2  or  AS4,  on  the  As  exposure  time,  and  on  the  temperature.  Fig. 
2.2.3  shows  the  As  3d  photoelectron  spectra  of  GaSb  layers  exposed  to  As  under  different 
conditions.  Under  the  present  flux  conditions,  for  an  exposure  time  of  10s  and  a  temperature  of 
400  °C,  no  appreciable  GaAs  is  observed.  The  formation  of  GaAs  on  the  surface  exhibits  a 
different  temperature  dependence,  depending  on  whether  As4  or  Ast  impinges  on  the  surface. 
In  particular,  when  As2  impinges  on  the  surface,  the  degree  of  exchange  increases  with 
temperature  (independent  of  exposure  time).  In  contrast,  when  As4  is  used,  GaAs  formation  is 
maximized  at  420°C.  In  addition,  exchange  increases  as  the  exposure  time  increases.  These 
different  temperature  dependences  imply  different  reaction  paths  for  anion  exchange 
depending  on  the  arsenic  species. 
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Figure  2.2.3:  As3d  photoelectron  peaks  acquired  for  GaSb  layers  exposed  to  As2  and  AS4 
for  different  times  and  at  different  temperatures. 


XPS  does  not  show  significant  Sb-for-As  exchange.  The  anion-to-Ga  atomic  percentage  ratios 
are  plotted  for  a  GaAs  layer  after  exposure  to  an  Sb2  flux  for  100s  at  550°C  (Fig. 2.2. 4a)  and  for 
GaSb  after  exposure  to  an  AS4  flux  for  10s  at  420°C  (Fig.  2.2.4b).  In  both  cases,  it  is  found  that 
the  total  anion-to-cation  ratio,  i.e.,  (Sb+As)/Ga,  is  higher  than  1,  indicating  an  excess  of  the 
anion  component  at  the  surface.  Indeed,  the  difference  between  the  two  systems  is  that  for  the 
Sb-for-As  anion  exchange  (Fig.  2.2.4a)  the  (Sb+As)/Ga  decreases  with  the  sampling  depth, 
indicating  that  the  anion  segregates  at  the  GaAs  surface,  while  for  the  As-for-Sb  anion 
exchange  (Fig.  2.2.4b),  the  (Sb+As)/Ga  ratio  increases  with  sampling  depth,  indicating  that  the 
anion  segregates  at  the  GaAsySbi.y/GaSb  interface.  Furthermore,  Fig.  2.2.4a  shows  that  it  is  the 
Sb/Ga  ratio  that  increases  at  the  outmost  surface,  providing  evidence  that  Sb  segregates  at  the 
GaAs  surface.  In  contrast,  Fig.  2.2.4b  shows  that  while  the  Sb/Ga  ratio  decreases  upon  anion 
exchange,  the  As/Ga  ratio  increases  with  sampling  depth  indicating  that  some  of  the 
incorporated  As  is  present  at  the  inner  interface  in  a  form  not  bonded  to  Ga.  Therefore,  for  the 
As  for  Sb  exchange,  the  trends  for  the  (As+Sb)/Ga  and  As/Ga  ratios  are  consistent  with  the 
segregation  at  the  inner  interface  of  AsSby  compounds. 
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Figure  2.2.4:  XPS  quantitative  ratios,  as  measured  from  the  peak  areas,  (a)  for  a  GaAs 
layer  as  grown  and  after  exposure  to  Sb2  for  100s  at  550°C,  and  (b)  for  a  GaSb  layer  as 
grown  and  after  exposure  to  AS4  for  30s  at  420°C. 


The  effects  of  temperature,  As  exposure  time,  and  As  species  on  the  extent  of  the  reactions  for 
As/GaSb  are  summarized  in  Figure  2.2.5.  It  can  be  seen  that  in  the  case  of  AS2  exposure,  no 
difference  is  found  between  the  10s  and  30s  exposure  times,  indicating  that  Ast  readily  reacts, 
and  exchange  is  not  under  kinetic/diffusion  control  under  the  present  conditions.  It  is  found 
that  AS4  is  less  reactive  than  AS2,  but  more  AsSby  is  formed.  However,  this  does  not  apply  for 
all  temperatures,  and  for  AS4  exposure,  a  critical  temperature  of  approximately  420°C  exists 
were  the  anion  exchange,  i.e.,  GaAs  fonnation,  is  favored  at  the  expense  of  the  AsSby. 
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Figure  2.2.5:  GaAs  relative  area  as  measured  from  the  fit  analysis  of  the  As3d  photoelectron  peaks  for 
GaSb  layers  exposed  to  As2  and  As4  for  10s  and  30s  at  various  temperatures. 

For  the  factors  investigated  in  this  study,  the  surface  temperature  has  the  most  significant 
impact  on  surface  morphology.  Figure  2.2.6  shows  AFM  topographic  images  of  GaAs  layers 
exposed  to  Sbi  for  100s  at  temperatures  of  420  °C  and  550  °C.  It  is  seen  that  surface 
modification  is  enhanced  as  the  temperature  is  increased.  The  Sb2  exposure  creates  an  oriented 
morphology,  similar  to  that  observed  for  AlGaAs  surfaces  exposed  to  indium  [24].  For  the  case 
of  As/Sb,  the  general  trend  is  that  surface  roughness  increases  with  temperature.  However,  a 
distinctly  different  morphology  is  observed  at  420  °C,  which  corresponds  to  the  temperature 
resulting  in  greater  GaAs  formation,  and  lower  AsSby  formation  as  previously  discussed.  In 
comparison  to  the  GaAs  exposure  to  Sb,  this  surface  shows  no  orientation  effect. 

Ellipsometric  analysis  of  two  samples  obtained  by  exposing  GaAs  surfaces  at  420  °C  and  550 
°C  to  Sb2  for  100s  are  indistinguishable  from  that  of  a  GaAs  substrate  with  a  native  oxide 
surface  layer.  This  indicates  that  significant  Sb  for  As  anion  exchange  does  not  occur.  In  order 
to  increase  the  modeling  accuracy,  spectra  were  acquired  at  various  angles  of  incidence  from 
55°  to  80°  for  the  550°C  sample,  and  were  fitted  simultaneously  to  a  superlattice  model 
GaAs/(GaSbyAsi-y/GaAs)2o/GaSb-oxide  where  fit  variables  were  the  layer  thickness’s  and  the 
exchanged  layer  composition. 
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Figure  2.2.6:  AFM  3D  topographical  images  for  GaAs  layers  exposed  to  Sb2  for  100s  at 
420°C  and  550°C,  and  for  GaSb  layers  exposed  to  As4  for  30s  at  various  temperatures. 
The  surface  roughness,  rms,  is  also  indicated.  Lateral  scale  is  lpmxlpm  for  every  image, 
but  5pmx5pm  for  last  image 


Figure  2.2.7  shows  the  best  fit  result  of  the  variable  angle  spectra.  These  results  indicate  that 
when  a  GaAs  layer  is  exposed  to  an  Sb  flux,  the  Sb  does  not  diffuse  into  the  GaAs  but,  instead, 
segregates  at  the  surface.  The  Sb  segregation  at  the  GaAs  surface  is  converted  into  GaSb 
and/or  GaSbyAsi.y  when  the  Ga  and  As  fluxes  are  initiated  for  the  overgrowth  of  the 
subsequent  GaAs  layer.  Furthermore,  no  diffusion  of  Sb  into  the  overgrown  GaAs  layer  is 
observed.  Thus,  an  abrupt  interface  is  achieved  in  (GaSbyAsi_y/GaAs)n  structures  formed  by 
letting  Sb2  impinge  on  a  GaAs  epitaxial  layer,  as  also  confirmed  by  the  intense  and  sharp 
interference  fringes  observed  by  XRD  analysis  of  (GaAs/GaSbyAsi_y)2o  structures  published 
previously. 


A  similar  approach  has  also  been  applied  to  the  characterization  of  the  As/Sb  anion  exchange 
and  interface  composition  in  (GaSb/GaAsySbi_y)2o  structures.  Figure  2.2.8  shows  a  comparison 
among  typical  spectra  of  the  real,  <si>,  and  imaginary,  <S2>,  parts  of  a  GaSb  substrate  and  the 
spectra  recorded  for  two  typical  (GaSb/GaAsySbi_y)2o  structures  obtained  by  exposure  of  the 
GaSb  epitaxial  layers  to  an  AS4  flux  for  10s  and  30s  at  a  temperature  of  420  °C.  The  spectra  of 
the  structures  differs  from  that  of  the  GaSb  substrate  with  native  oxide,  mainly  in  the  region  of 
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the  Ei  and  E1+A1  interband  critical  points  (CPs)  that  are  strongly  sensitive  to  material 
composition.  In  particular,  for  short  As  exposure  time,  an  increase  of  the  Ei  and  E1+A1  CPs 
amplitude  is  observed,  while  it  decreases  with  respect  to  that  of  the  GaSb  substrate  with  native 
oxide  for  longer  As  exposure  time.  This  indicates  that  the  exchange  occurs  to  a  different  extent 
depending  on  exposure  time.  A  simple  model  assuming  the  formation  of  a  ternary  alloy  does 
not  fit  the  amplitude  of  the  Ei  and  E1+A1  CPs. 


Figure  2.2.7:  SE  spectra  acquired  at  various  angles  of  incidence  from  55°  to  80°  for  a  typical  20-period 
(GaSbyAsi_y)/GaAs  obtained  by  GaAs  layer  exposure  to  an  Sb2  flux  for  100s  at  550  °C. 


A  good  fit  is  obtained  by  a  GaSb  buffer  layer  thickness  of  about  270 A,  and  an  exchanged  layer 
with  a  total  thickness  of  approximately  19  A  that  includes  a  very  thin  interface  layer  containing 
some  AsSb.  The  contribution  of  AsSb  may  be  responsible  for  the  increase  of  the  Ei  and  Ei+Ai 
CPs  amplitude  found  for  the  sample  exposed  for  10s  to  arsenic.  The  dependence  of  the 
thickness  of  the  exchanged  layer  and  of  the  AsSby  segregation  on  the  As  species,  temperature, 
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and  As  exposure  time  is  summarized  in  Fig.  2.2.9.  It  is  seen  that  the  thickness  of  the  exchanged 
layer  is  15-16A±lA  for  an  As  exposure  time  of  10s  and  is  16-18A±lA  for  an  AS2  exposure 
time  of  30s.  These  data  agree  with  the  XPS  data  in  Fig.  2.2.6,  where  a  difference  of  1-2  A 
cannot  be  resolved.  Furthennore,  the  thickness  of  the  exchanged  layer  increases  with 
temperature  for  As2,  while  it  has  a  maximum  at  420  °C  for  AS4.  The  measured  thickness 
dependence  can  be  related  to  AsSby  formation  and  segregation  at  the  GaAs/GaSb  interface 
reported  in  Fig.  2.2.9b.  In  fact,  the  lower  the  AsSby  segregation,  the  greater  the  thickness  of  the 
exchanged  layer,  indicating  that  the  AsSby  segregation  inhibits  anion  exchange. 
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Figure  2.2.8:  Typical  spectra  of  the  real,  <ei>,  and  imaginary,  <82>,  parts  of  a  GaSb 
substrate  and  the  SE  spectra  recorded  for  two  (GaSb/GaAsySbi-y)2o  SLs  obtained  by 
exposure  of  the  GaSb  epitaxial  layers  to  an  As4  flux  for  10s  and  30s  at  the  temperature  of 
420°C. 
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Figure  2.2.9:  Temperature  and  time  dependence  of  (a)  the  total  exchanged  layer  and  (b)  AsSby  interface 
segregation  as  derived  by  SE  analysis  of  (GaSb/GaAsySbl  y)20  SLs  obtained  by  As2  and  As4  soak  of  GaSb 
layers. 


4.2.3.  P/As  Reactions 

The  formation  of  Ga-P  bonds,  indicative  of  P-for-As  exchange,  is  observed  in  the  P  2p 
photoelectron  spectra  acquired  for  GaAs  surfaces  exposed  to  a  P2  flux  at  different  temperatures 
and  exposure  times.  Figure  2.2.10  shows  the  P/Ga  ratio  as  calculated  from  the  Ga  3d  and  P  2p 
peak  areas  for  different  temperatures  and  P  exposure  times  It  is  found  that  the  higher  the 
temperature  and  exposure  time,  the  larger  the  P  surface  coverage,  which,  however,  is  lower  for 
the  As-rich  GaAs  surface  than  for  the  As-lean  surface.  This  indicates  that  more  As  on  the  GaAs 
surface  hinders  P  surface  coverage  and  exchange.  At  high  temperature  (520  °C)  and  long  P 
exposure  time,  the  two  surfaces  behave  similarly  because  of  the  thermal  desorption  of  the 
excess  of  As  that  acts  as  a  surface  barrier  for  the  P  uptake.  The  As-rich  reconstruction  of  the 
GaAs  surface  reduces  P  incorporation,  in  agreement  with  previous  work  [25], 
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Figure  2.2.10:  Surface  P/Ga  ratio  derived  from  area  ratio  of  the  P2p  and  Ga3d  XPS  peaks  for  As-rich  and 
As-lean  GaAg  layers  exposed  to  P2  at  various  temperatures  for  10s  and  30s.  The  XPS  take-off  angle  is  20°. 
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Figure  2.2.11:  Spectra  of  the  imaginary,  <e2>,  part  of  the  pseudodielectric  functions  acquired  for 
(GaAs/GaPyAsi_y)10  SLs  by  exposure  of  GaAs  to  a  P  flux  at  420°  for  10s  and  50s  and  at  520°C  for  50s.  Inset: 
Detail  of  the  Ej  and  E^AE!  CPs.  The  corresponding  best-lit  models  at  bottom 


Figure  2.2.11  shows  typical  ellipsometric  spectra  of  the  imaginary,  <S2>,  part  of  the 
pseudodielectric  function  of  (GaAs/GaPyAsi_y)io  structures  obtained  by  exposure  of  GaAs  to  a 
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P  flux  for  different  times  and  at  different  temperatures.  The  corresponding  best  fit  BEMA 
models  are  also  shown,  as  well  as  details  of  the  Eo,  Ei  and  E1+AE1  interband  critical  CPs.  The 
details  of  Ei  and  E1+AE1  CPs  (see  inset)  shows  a  blue-shift  of  the  CPs  and  a  gradual  reduction 
with  the  increase  of  surface  temperature  and  phosphorous  exposure  time.  The  dielectric 
function  of  GaPyAsi.y  ternary  alloys  show  a  gradual  evolution  of  the  characteristic  interband 
critical  points  with  increasing  phosphorous  content  from  that  of  GaAs  to  that  of  GaP  [26],  We 
derived  the  composition  of  the  exchanged  GaPyAsi_y  layer,  which  ranges  from  y  =  0.1  -  0.2  for 
the  investigated  temperature  and  phosphorous  exposure  time  conditions.  It  is  found  that  the 
GaPyAsi-y  stoichiometry  mainly  depends  on  the  GaAs  surface  reconstruction  and  on  the 
temperature.  In  particular,  for  the  As-rich  GaAs  surface  reconstruction,  a  composition  of  0.1  is 
found.  In  contrast,  for  the  As-lean  GaAs  surface  reconstruction,  the  phosphorous  composition 
increases  to  0.2.  The  thickness  of  the  exchanged  layer  increases  with  phosphorus  exposure  time 
and  temperature.  In  particular,  the  modeling  analysis  shows  that  for  As-rich  GaAs  surface 
reconstruction  and  for  T  <  520  °C,  P-for-As  anion  exchange  occurs  over  approximately  100 A 
yielding  a  graded  layer  from  the  ternary  GaP0.iAso.9  at  the  surface  to  GaAs.  A  GaAs  layer 
approximately  100 A  thick  that  does  experience  exchange  is  still  present.  In  contrast,  for  the 
highest  temperature  of  520  °C  and  especially  for  the  As-lean  GaAs  samples,  P  for  As  exchange 
and  intermixing  occurs  over  the  entire  200A  GaAs  epitaxial  layer.  The  data  for  the  exchanged 
layer  are  in  good  agreement  with  XRD  analysis. 


Figure  2.2.12:  Dependence  of  integrated  P  content  on  P  soak  time,  as  estimated  by  analysis  for  As-rich  (PAs 
=  4><10"6  Torr)  and  As-lean  (PAs  =  4xl0'6  Torr)  GaAs  epitaxial  layers  exposed  to  the  P2  flue  at  various 
temperatures. 

The  integrated  phosphorus  content  has  been  calculated  and  reported  in  Fig.  2.2.12  as  a  function 
of  phosphorus  exposure  time  for  various  temperatures.  It  is  found  that  the  phosphorus  content 
is  higher  for  the  As-lean  GaAs  epitaxial  layer  than  for  the  As-rich  GaAs,  and  in  both  cases 
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slightly  increases  with  time  ranging  from  10s  to  50s.  The  phosphorus  content  quickly  increases 
in  the  first  few  seconds  of  exposure  to  phosphorus,  indicating  a  very  fast  in-diffusion  of  P  into 
GaAs.  Therefore,  anion  exchange  does  not  occur  in  the  diffusion-controlled  kinetic  regime 
under  the  present  conditions. 

Figure  2.2.13  shows  typical  XRD  rocking  curves  obtained  from  structures  exposed  to  P2  for  30 
seconds  and  varying  temp  for:  a)  As  rich  surfaces,  and  b)  As-lean  surfaces.  For  lower 
temperature  and  exposure  time  and  As-rich  surfaces,  the  x-ray  diffraction  data  indicate  the 
formation  of  a  well-ordered  layers  with  many  positive  and  negative  superlattice  peaks.  The  P 
composition  is  calculated  from  the  separation  of  the  0th  order  reflection  and  the  substrate  peak. 
The  calculations  assume  that  ah  of  the  strain  is  elastically  accommodated  -  an  assumption  that 
is  not  equally  accurate  for  ah  the  samples,  but  provides  an  observation  of  trends  in  the  data  for 
comparison  to  SE  and  XPS.  The  dramatic  difference  in  the  phosphorus  incorporation  is  clear  in 
the  comparison  of  the  x-ray  for  the  As-rich  and  lean  surfaces.  The  absence  of  well-ordered 
negative  peaks  in  Fig.  2.2.13b  indicates  the  much  greater  incorporation  of  P  in  these  structures 
and  the  greater  degree  of  material  strain  and  disordering.  A  comparison  of  x-ray  data  and  the 
simulations  used  to  fit  the  profiles  are  shown  in  Fig.  2.2.14.  The  profile  simulated  for  SE  and 
XRD  are  shown  for  comparison  and  indicate  the  high  consistency  obtained  for  the  two 
different  measurements  and  simulations.  Finally,  the  predictions  of  a  neural  network  process 
model  that  was  derived  from  the  XRD  data  are  shown  in  Fig.  2.2.15.  The  surfaces  show  the 
dependence  of  P  content  on  time  and  temperature  and  also  show  the  measured  variant  of  the 
(2x4)  surface  reconstruction  for  the  ranges  of  growth  conditions. 

Surface  topography  was  measured  by  atomic  force  microscopy  (AFM).  Fig.  2.2.16  shows  the 
range  of  morphologies  observed  for  the  samples  synthesized  at:  a)  low  temperature  and 
exposure  time  and  As-rich  surface,  and  b)  high  temperature  long  exposure  time  and  As-lean 
surface.  At  lower  temperature,  the  morphology  is  characterized  by  the  appearance  of  mounds 
which  are  typical  of  low  temperature  growth  of  GaAs.  At  high  temperature  the  surface  exhibits 
a  cross-hatch  morphology  indicative  of  misfit  dislocation  formation  and  some  degree  of 
relaxation  of  the  strain  realized  in  the  GaAsP  layer.  It  is  not  clear  if  the  strain  and  relaxation  is 
uniform  throughout  the  depth  of  the  layer  or  gradually  increases  with  the  superlattice  period. 
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(a) 


Figure  2.2.13:  XRD  data  (004)  for  30  sec  P  exposures  for  As-rich  surface  (a)  and  As-lean 

surfaces  (b) 
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Figure  2.2.14:  XRD  data  and  fit  with  simulation  for  XRD  (a)  and  SE  (b)  showing 

simulated  structure  comparison 
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Figure  2.2.15:  Neural  network  process  models  of  P  composition  (from  XRD  analysis)  as  a 
function  of  temperature  and  exposure  time  for  (a)  As-lean  and  (b)  As-rich  surfaces. 
Variants  of  the  (2x4)  reconstruction  are  also  shown 
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Figure  2.2.16:  Topographic  AFM  iunages  for  P/GaAs  (a)  low  temperature  and  exposure 
time  and  As-rich  surface  and  (b)  high  temperature  long  exposure  time  and  As-lean 
surface 


4.2.4.  Chemical  Models  for  Sb/As,  As/Sb,  and  P/As  Surface  Reactions 

The  Sb/GaAs  is  the  most  simple  system  examined  and  is  characterized  by  Sb  adsorption  and 
subsequent  incorporation.  We  therefore  model  it  as: 

£  £)  , 

Sb(g)  < — »  Sb(ads)  (2)  ►Sb(s)  +  GaAs-5^  no  reaction  (17) 

Very  low  tV 

where  (1)  is  the  surface  chemisorption  of  antimony,  (2)  is  the  diffusion  step  and  (3)  is  the 
reactive  step  with  the  corresponding  rate  constants. 

Surface  reactivity  requires  overcoming  the  energetic  barrier  for  surface  to  bulk  atomic  in- 
diffusion  ( E  s^b )  and  for  bulk  to  surface  exchanged  anion  out-diffusion  (Eb^s).  For  As  on 

GaSb,  the  driving  force  Es^h  -  Ebs^s  >  0  is  because  of  the  faster  in-diffusion  of  As;  however, 

the  very  low  out-diffusion  of  Sb  causes  AsSby  to  segregate  at  the  GaAs/GaSb  interface. 
Therefore,  the  in-diffusion  of  As  is  not  rate  limiting,  and  as  demonstrated  by  the  data  presented 
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for  AS2  which  indicates  that  the  exchange  is  not  kinetically  limited,  we  conclude  that  a  simple 
diffusion-based  model  cannot  explain  the  different  products  and  reactivity  observed  for  AS2 
and  AS4.  For  this  reaction,  the  following  steps  should  be  considered:  (1)  chemisorption  of  the 
As  on  the  GaSb  surface;  (2)  in-diffusion  of  As  through  the  GaSb  layer  and  through  the  product 
GaAsySbi.y  layer;  (3-6)  chemical  reactions  at  the  GaAs/GaAsySbi_y  surface  involving  As  for  Sb 
exchange  with  the  formation  of  GaAs  and  of  the  byproduct  involving  the  exchanged  Sb,  i.e., 
AsSby;  (7)  outward  diffusion  of  the  Sb  products  and  desorption  or,  as  in  the  present  case, 
interface  segregation.  Based  on  these  assertions  and  the  data,  the  chemical  model  for  the 
interaction  of  As2  and  As4  with  GaSb  can  be  described  by: 


Favoured  by  As2 
k \ 

As  (g)  As  (ads)-—*  As  (s)  +  GaSb 

(1)  (2) 


Favoured  by  As4 


GaAs„Sb  +  Sb  Js) 


+As 

GaSbj.y 
Isoelectronic  AsSb  formation 


(18) 


where  Kacis  is  the  chemsiorption  constant,  D  is  the  diffusion  coefficient  for  the  indiffusion  of  As 
atoms,  keXch  and  kiso  are  the  true  rate  constants  of  the  reactive  step  at  the  GaSb  surface  yielding 
anion  exchange  and  isoelectronic  formation,  respectively,  and  kr  and  kr’  are  the  rate  constants 
for  the  recombination  of  the  exchanged  Sb  and  of  the  As  vacancy  on  a  Ga  site  (GaSbi_y), 
respectively.  According  to  the  proposed  chemical  model,  the  different  temperature 
dependences  observed  for  the  reaction  for  As4  and  As2  species  implies  different  reaction 
pathways.  For  both  species,  the  final  products  are  GaAsySbi.y  and  AsSby  that  form  in  steps  (3) 
+  (6)  and  (4)  +  (5)  in  the  case  of  As2  and  As4,  respectively.  Furthermore,  these  reaction 
pathways  offer  an  explanation  for  the  different  concentrations  of  the  reaction  products  of  the 
exchange  and  AsSby  formation  for  the  As2  and  As4  species.  In  particular,  for  the  As2,  the  larger 
GaAsySbi.y  formation  results  from  the  first  order  kinetics  being  kexch  >  kr.  For  the  As4  species, 
the  observed  temperature  dependence  through  a  maximum  for  the  GaAsySbi.y  formation  could 
be  explained  by  applying  the  formalism  of  the  topochemical  reactions  for  a  second  order 
kinetics,  i.e. 


Sb  I 

- =  Ko  •  K  ■  el  •  (19) 

at 

where  ©as  is  the  coverage  of  As  at  the  reactive  interface  and  SGaSb  is  the  GaSb  surface  area 
available  for  the  anion  exchange.  The  above  equation  takes  into  account  the  different 
temperature  dependence  of  the  k,s0  and  kf  rate  constants,  and  of  the  ©as  as  well. 
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Previous  studies  [27]  have  reported  that  As  for  Sb  exchange  is  self-limiting  because  after  a  few 
monolayers  of  GaAs  formation,  further  exchange  is  inhibited.  According  to  the  proposed 
model,  the  explanation  of  the  self-limiting  nature  is  the  interface  segregation  of  AsSby  that  acts 
as  a  sink  for  As,  and  reduces  SoaSb-  This  competition  between  AsSby  formation  and  anion 
exchange  may  explain  reports  of  disordered,  defective,  and  phase  separated  GaAs/GaSb 
superlattices  using  AS4,  while  structures  grown  with  AS2  show  improved  crystalline  quality. 

The  case  of  P/As  exchange  is  examined  in  the  context  of  surface  reconstruction  and  associated 
As  population.  The  experimental  results  presented  above  indicate  that  an  increased  surface 
population  of  As  acts  as  a  barrier  to  the  formation  of  GaPyAsi.y.  Therefore,  P-for-As  exchange 
is  best  described  in  terms  of  the  concomitant  and  competitive  reactions  of  both  As  and  Ga  with 
incident  P.  The  rates  of  these  two  reactions  determine  the  rate  of  the  overall  process  and  the 
composition  of  the  layer  created  by  phosphorus  in-diffusion.  In  addition,  the  reactions  depend 
on  the  temperature.  In  order  to  clarify  these  effects,  we  begin  from  a  consideration  of  the 
different  factors  controlling  the  reactivity  of  the  P  and  As.  From  a  thermodynamic  point  of 
view,  AsP  (4.49  eV)  and  GaP  bond  strengths  (2.38  eV)  are  higher  than  GaAs  (2.17  eV).  Hence, 
the  higher  chemical  stability  of  products  is  the  driving  force  for  the  reaction.  In  the 
temperature  range  investigated,  the  desorption  rate  constant  of  arsenic  is  much  smaller  than 
that  of  phosphorus.  The  activation  energies  for  arsenic  and  phosphorus  desorption  from  GaAs 
and  GaP  are  2.5  eV  and  1.86  eV,  respectively  [27].  The  diffusion  coefficient  of  P  into  GaAs 
( 1  O'7- 1 0  s  cmV1)  is  higher  than  the  self-diffusion  coefficient  of  As  in  GaAs  [28].  This  means 
that  the  P  in-diffusion  is  very  fast  and  not  rate  limiting.  Given,  the  lower  coefficient  of  self- 
diffusion  of  As  in  GaAs,  through  vacancies  that  are  saturated  with  in-diffusing  phosphorus,  the 
out-diffusion  of  As  could  be  a  slow  process.  However,  we  do  not  observe  by  XPS  or  SE  any 
interface  segregation  of  As  resulting  from  the  P-for-As  anion  exchange.  From  this,  it  can  be 
inferred  that  P  can  easily  react  with  As  either  in  the  solid  or  on  the  surface  at  temperatures 
yielding  PAs,  which  has  a  high  vapour  pressure.  On  the  basis  of  the  earlier  considerations,  we 
can  begin  to  outline  a  chemical  model  for  P-for-As  exchange  at  GaAs  surfaces.  The  model 
should  include  the  supply  of  P  to  the  GaAs  surface  where  P  adsorbs  and  the  chemical  barrier 
effects  of  the  surface  excess  of  As.  We  will  assume  a  surface  reactive  step  involving  the  As 
and  the  P,  yielding  volatile  AsP  that  getters  As  from  the  GaAs  surface  and  makes  Ga  sites 
available  for  exchange.  Due  to  the  fast  P  in-diffusion  into  GaAs,  P  for  As  exchange  proceeds 
through  the  bulk  to  yield  GaP  bond  formation,  and  As  atom  can  slowly  diffuse  outward  and/or 
react  quickly  with  P  yielding  PAs.  All  the  above  steps  can  be  included  in  the  following 
chemical  model  for  the  P  for  As  interactions: 
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Favoured  by  As-lean 
surface  reconstruction 


K, 


D 


P  (g)  «=►  P  (ads)  — ^>P  (s)  +  GaAs 


(i) 


Favoured  by  As-  rich 
surface  reconstruction 


Anion  Exchange 


GaPyAsj.y  +  yAs  (s) 


(6) 


+P 


GaAsj 


+  yAsPt 


Isoelectronic  AsP  formation 

(at  GaAs  surface) 


(20) 


where  Kac/S  is  the  chemsiorption  constant  for  P  atoms,  D  is  the  diffusion  coefficient  for  the  in- 
diffusion  of  P  atoms,  kexch  and  klS0  are  the  true  rate  constants  of  the  reactive  step  at  the  GaAs 
surface  yielding  anion  exchange  and  an  As-vacancy  by  isoelectronic  formation,  respectively 
and  k’  is  the  rate  constant  for  the  reaction  of  P  with  the  As  vacancy  on  a  Ga  site  (GaAsi_y) 
yielding  the  ternary  GaPyAsi.y  compound.  According  to  the  above  chemical  model,  the 
formation  of  GaPyAsi.y,  occurs  through  a  first  order  kinetics  on  the  As-lean  GaAs,  while  it  is  a 
two  step  reaction  on  the  As-rich  GaAs.  Thus,  the  above  experimental  data  have  shown  that, 
independent  of  As  surface  reconstruction  of  GaAs,  the  P  for  As  exchange  process  is  not 
controlled  by  the  P  diffusion.  This  indicates  that  the  reactive  step  (3-4)  takes  place  in  the 
kinetic  regime,  i.e.,  the  rate  of  the  reactive  step  depends  on  the  GaAs  sites  available  for  P 
adsorption.  Hence,  an  excess  of  As  inhibits  the  formation  of  GaPyAsi_y  by  scavenging  P 
through  the  formation  of  PAs. 


2.3.  Hybrid  Neural  Network  Modeling  of  Anion  Exchange 

Using  some  simplifying  assumptions,  a  hybrid  neural  network  model  was  constructed  by 
characterizing  the  growth  of  GaAsi_yPy/GaAs  strained  layer  superlattices  (SLS)  grown  on  (001) 
GaAs  substrates  by  molecular  beam  epitaxy.  These  heterostructures  were  formed  by  the  P? 
exposure  of  an  As-stabilized  GaAs  surface,  and  ex-situ  high  resolution  x-ray  diffraction 
(HRXRD)  was  performed  to  determine  the  phosphorus  composition  at  the  interfaces.  A  first- 
order  kinetic  model  was  then  developed  to  describe  the  mechanisms  of  anion  exchange,  surface 
desorption,  and  diffusion.  A  semi-empirical  hybrid  neural  network  was  used  to  estimate  the 
parameters  of  the  kinetic  model  and  analyze  the  microscopic  processes  occurring  at  the 
interfaces  of  the  mixed  anion  III-V  heterostructures.  The  phosphorus  diffusion  process  in  GaAs 
was  estimated  to  have  a  diffusion  coefficient  of  D  =  7x  10'15  exp  (-0. 1 1  /kBTs)  curs' 1  for  samples 
with  Pas4  =  4x  1 O'6  torr  and  exhibited  enhanced  phosphorus  intermixing  for  samples  with  lower 
As-stabilizing  fluxes. 

The  investigation  proceeded  by  conducting  a  D-optimal  experiment  designed  to  systematically 
characterize  the  growth  of  10-period  GaAsi_yPy/GaAs  strained  layer  superlattice  (SLS) 
structures.  The  incorporation  of  phosphorus  in  the  SLS  structures  was  achieved  by  the  P2 
exposure  of  an  As-stabilized  GaAs  surface  and  an  immediate  overgrowth  of  GaAs.  The  D- 
optimal  experiment  implemented  for  this  investigation  provided  a  basis  for  modeling  the 
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phosphorus  composition  at  the  interfaces  of  the  SLS  structures  as  a  function  of  substrate 
temperature  (Ts),  phosphorus  exposure  time  (4xP),  and  As-stabilizing  flux.  The  growth 
conditions  and  factor  settings  for  the  D-optimal  experiment  are  provided  in  the  Table  2.3.1 
below. 


Factor  Description 

Symbol 

Settings 

Substrate  Temperature 

Ts 

420  to  520°C 

P2  Exposure  Time 

texp 

10  to  50s 

P2  Exposure  Flux  (BEP) 

p? 

r2 

4.5x1  O'6  Ton- 

AS4  Stabilizing  Flux  (BEP) 

Pas4 

2xl0’6  Torr,  4xl0‘(1  Ton 

Growth  was  performed  on  undoped  (001)  GaAs  epi-ready  substrates  in  a  Riber  2300  solid- 
source  molecular  beam  epitaxy  system  equipped  with  a  valved  As  cracker  and  a  three-zone 
valved  P  cracker.  The  surface  reconstructions  at  the  interfaces  of  the  SLS  structures  were 
observed  during  P2  exposure  of  GaAs  using  in-situ  reflection  high-energy  electron  diffraction 
(RHEED).  The  phosphorus  composition  present  at  the  interfaces  of  the  GaAsi.yPy/GaAs  SLS 
structures  was  determined  by  performing  standard  HRXRD  #/2i9  lincscans  around  the  substrate 
(004)  reflections. 


The  following  kinetic  model  was  derived  to  quantitatively  describe  the  microscopic  processes 
that  occur  during  the  complete  growth  sequence  for  the  SLS  structures. 


y(Ts .  Lxp ,  o)  =  sOt [l  -  exp(-  texp  /  r)]  x  [l  +  2 ^Dtexp  /  n 


(21) 


A  semi-empirical  hybrid  neural  network  was  subsequently  developed  that  incorporates  this 
first-principles  kinetic  model  with  a  multilayer  feed-forward  perceptron  neural  network.  The 
hybrid  neural  network  designed  to  predict  anion  intermixing  for  the  GaAsP/GaAs 
heterostructures  is  illustrated  below.  The  neural  network  component  of  the  hybrid  model  has 
the  MBE  process  conditions  as  its  inputs.  The  outputs  of  the  neural  network  component  are  the 
unknown  parameters  required  to  implement  the  kinetic  model. 
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Growth  Kinetic  Predicted 


Composition 


Figure  2.3.1:  Hybrid  neural  network  used  to  estimate  the  kinetic  parameters  for  growth  of  the  As/P 
heterostructures. 


Evaluation  of  the  hybrid  neural  network  model  was  performed  in  terms  of  the  root  mean 
squared  error  (RMSE).  The  hybrid  neural  network  implemented  for  samples  with  an  As- 
stabilizing  flux  P As4  =  4x1  O'6  Torr  demonstrated  a  training  RMSE  of  0.028%  and  a  prediction 
RMSE  of  0.574%.  The  figure  below  illustrates  the  measured  phosphorus  composition  versus 
the  predicted  composition  for  these  samples,  indicating  the  accuracy  of  the  hybrid  neural 
process  model. 


Figure  2.3.2:  Experimental  phosphorus  composition  versus  neural  network  predictions  for  As4  =  4x10  6 
torr. 
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Similar  results  were  obtained  for  the  samples  with  an  As-stabilizing  flux  Pas4  =  2xl0'6  torr. 
The  measured  phosphorus  composition  versus  the  predicted  composition  for  these  samples  is 
illustrated  in  the  figure  below.  The  kinetic  parameters  derived  are  provided  in  Table  2.3.2. 
Using  these  estimated  parameters,  the  hybrid  model  can  predict  the  phosphorus  composition  of 
the  mixed  arsenic/phosphorus  III-V  heterostructures  accurately. 


Figure  2.3.3:  Experimental  phosphorus  composition  versus  neural  network  predictions  for  As4  =  2x10 6 
torr. 


The  hybrid  model  is  important  for  analyzing  the  microscopic  processes  occurring  at  the 
interfaces  of  the  GaAsP/GaAs  heterostructures.  Several  fundamental  relationships  are  evident 
upon  comparing  the  kinetic  parameters  for  the  two  distinct  As-stabilizing  fluxes.  The  higher  As 
flux  corresponds  to  atomic  surface  structures  with  increased  surface  arsenic  coverage.  It  is 
apparent  from  the  model  predictions  that  the  phosphorus  sticking  coefficient  for  samples  with 
Tas4  =  4x1 0‘6  torr  is  smaller  than  those  with  Pas4  =  2x1  O’6  torr.  This  is  due  to  the  relative 
increase  in  As  surface  coverage  corresponding  to  stable  surface  structures  such  as  the  c(4x4) 
reconstruction.  Experimental  estimates  of  group  V  sticking  coefficients  predict  values  ranging 
from  0.25  to  0.75  cm2.  The  sticking  coefficient  of  P2  on  the  (001)  GaAs  surface  is  at  a 
maximum  when  the  group-III  is  present  on  the  surface.  Therefore,  the  phosphorus  sticking 
coefficient  values  of  -0.36  cm2  obtained  correspond  well  to  the  non-growing  (no  incident 
group-III)  anion  exposure  conditions. 
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Table  2.3.2:  Kinetic  parameters  from  neural  network  model. 


Kinetic  Parameters  - 

-  AS4  =  4xl0'6  Torr 

s 

Ad 

Ed 

Do 

Ea 

(s) 

(eV) 

(mol/s) 

(eV) 

0.358 

6.412 

2.985 

17.358 

7xl0'15 

0.11 

Kinetic  Parameters  - 

-  AS4  =  2xl0'6  Torr 

s 

Ad 

Ed 

Do 

Ea 

wmm 

(s) 

(eV) 

(eV) 

0.369 

3.665 

2.427 

26.994 

lxlO"14 

0.05 

For  phosphorus  surface  desorption,  a  lower  activation  energy  and  increased  reaction  frequency 
(shorter  surface  lifetime)  were  observed  for  the  PA s4  =  2x1  O'6  torr  samples.  These  results 
indicate  an  increase  in  phosphorus  desorption  for  the  Pas  =  2xl0"6  torr  samples.  It  should  be 
noted  that  a  decrease  in  As  pressure  results  in  a  less  As-rich  surface,  and  group  V  desorption 
begins  to  occur  at  lower  substrate  temperatures.  Therefore,  the  a  phase  of  the  GaAs  (001)- 
(2x4)  surface  structure,  obtained  at  Ts  =  520  °C  for  the  Pas4  =  2x1  O'6  torr  samples,  has 

increased  group- V  desorption  compared  to  the  /T2x4)  surface  structure,  present  at  Ts  =  520  °C 
for  the  PAs4  =  4x1  O'6  torr  samples.  These  results  support  the  theory  that  surface  desorption  is 
governed  by  the  species  of  atoms  and  their  orientation  (i.e.,  reconstruction)  on  the  surface  and 
contributes  to  the  characteristic  nonlinear  behavior  of  group  V  incorporation. 


2.4.  InAs-AlSb  HEMT  Structures:  Impact  of  Inverted  Interface  Formation 

The  InAs-AlSb  high  electron  mobility  transistor  (HEMT)  is  of  critical  importance  for  future 
millimeter  wave  and  high  speed  electronic  devices.  We  have  achieved  record  300K  electron 
mobility  in  InAs  channel  structures.  We  have  applied  various  growth  techniques  to  the  critical 
mixed-anion,  As/Sb,  inverted  interface  to  discover  the  impact  on  device-relevant  properties. 

Experiments  were  perfonned  so  that  the  effect  of  growth  conditions  on  the  interface  formation 
could  be  examined.  The  statistically  designed  experiment  conditions  and  results  are  provided 
in  Table  2.4.1.  Disregarding  the  outlier,  these  experiments  produced  test  structures  with 
average  300  °K  mobility  above  20,000  cm/V/s.  They  also  produced  low  temperature  mobility 
that  encompassed  a  range  of  90,000  cm2/V/s. 
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Table  2.4.1:  300  °K  and  77  °K  Mobility  Results 


INDIUM 
BARRIER  (A) 

CHANNEL  TEMP 
=  460  °C+ 

77  °K  MOBILITY 

(cm2/V/s) 

300  °K  MOBILITY 

(cm2/V/s) 

15 

-15 

95038 

24773 

3 

30 

106792 

29356 

17.48 

7.5 

135367 

24259 

15 

30 

80257 

25716 

3 

-15 

91439 

28520 

9 

39.31 

98921 

24797 

1.3 

7.5 

101738 

31954 

9 

7.5 

137729 

28188 

Neural  network  process  models  were  developed  for  each  response.  Growth  conditions 
(temperature,  indium  barrier  thickness)  were  inputs  to  the  neural  network,  and  the  responses 
(electron  mobility  at  300  °K  and  77  °K)  were  outputs.  The  network  structures  and  results  are 
provided  in  Table  2.4.2.  Plots  of  the  growth  conditions  and  responses  generated  by  the  process 
models  are  provided  in  Figure  2.4.1. 


Table  2.4.2:  Growth  condition  process  models 


GC  300 

GC  77 

°K 

°K 

NN  Structure 

2-2-1 

2-2-4- 1 

Training  error 
(%) 

0.23 

0.48 

Testing  error 
(%) 

5.35 

6.55 
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Figure  2.4.1:  Models  based  on  growth  conditions  (a)  300  °K  mobility  (b)  77  °K  mobility 


It  is  well  known  that  In-Sb  bond  formation  must  be  forced  at  the  inverted  interface  of  InAs- 
AlSb  in  order  to  achieve  high  mobility  in  the  InAs  channel.  The  reasons  for  this  are  still 
somewhat  unclear,  although  clues  suggest  that  exchange  and  segregation  play  a  key  role. 
Therefore  we  choose  to  investigate  the  conditions  of  substrate  temperature,  which  is  intimately 
related  to  segregation  and  exchange,  and  In  deposition  on  the  Sb-saturated  AlSb  surface.  The 
room  temperature  mobility  model  derived  from  the  experimental  data  indicates  a  strong 
correlation  with  substrate  temperature.  The  electron  mobility  is  inversely  proportional  to 
substrate  temperature  during  formation  of  the  interface.  Conversely,  the  deposited  indium 
thickness  does  not  appear  to  have  an  effect  on  the  room  temperature  mobility.  Depositing 
multiple  monolayers  of  indium  at  the  interface  seems  illogical  since  it  leaves  the  interface 
metal  rich,  except  the  results  do  not  invalidate  these  experiments.  It  may  be  that  room 
temperature  electron  transport  measurements  are  not  sensitive  enough  to  explore  the  impact  of 
both  indium  and  temperature  on  the  interface  formation.  However,  it  is  certain  that  substrate 
temperature  plays  a  role  in  the  amount  of  indium  retained  on  the  growth  surface,  ie.  some  of 
the  InSb  may  well  be  desorbed  from  the  surface  thereby  modifying  the  Sb/As  profile  at  the 
interface. 

Low-temperature  electron  transport  is  sensitive  to  the  quality  of  the  interface  [29].  Interfacial 
roughness  scattering  has  been  responsible  for  the  decrease  in  electron  mobility  at  low 
temperature.  The  low-temperature  electron  transport  presented  demonstrate  a  correlation  with 
both  temperature  and  indium  deposition.  The  low-temperature  model  (Fig.  26b)  reveals 
maximums  in  two  particular  regions:  (1)  more  indium  deposited  at  higher  temperature,  and  (2) 
thinner  indium  layer  at  lower  temperature.  It  is  possible  that  there  is  little  difference  in  the 
quality  of  the  interface  between  these  two  regions.  In  region  (1),  indium  desorbs  at  higher 
temperatures,  thus  creating  an  interface  similar  to  those  formed  in  region  (2).  However,  when 
the  indium  barrier  is  thick  at  lower  temperature,  multiple  layers  of  indium  can  exist,  which  are 
converted  to  InAs  by  the  ensuing  As  flux.  The  model  also  illustrates  that  when  the  indium 
barrier  is  thin,  low  temperature  mobility  is  inversely  proportional  to  increasing  temperature, 
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suggesting  that  interface  is  affected,  possibly  roughened,  as  indium  desorbs  from  the  surface. 
When  indium  is  desorbed  so  as  not  to  protect  the  underlying  Sb  layer,  anion  exchange  can 
occur,  thus  creating  AlAs-bonds  at  the  interface.  The  low-temperature  model  is  effective  in 
illuminating  indium’s  role  as  a  function  of  substrate  temperature  on  interface  fonnation. 
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